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PREFACE 


Large  and  small  calibre  guns  have  long  been  important  tools  of 
warfare  and  still  are.  A  detailed  understanding  of  the  gas  dynamics 
in  a  discharged  gun  could  be  of  great  help  in  the  design  of  new  guns, 
shells  or  charges.  Such  knowledge  is  also  important  to  efforts  to 
suppress  or  detect  gun  smoke,  gun  flash,  and  gun  blast.  The  work  re¬ 
ported  on  in  this  memorandum  was  prompted  by  a  need  for  quantitative 
data  on  temperatures  and  densities  in  the  gun  exhaust  pertinent  to  de¬ 
tection  by  radar  or  ■'nfrared  scanning  devices.  Dr.  William  G.  McMillan 
first  directed  our  attention  to  this  problem,  and  we  are  further  in¬ 
debted  to  him  for  suggesting  sources  for  experimental  data. 

The  methods  are  developed  here  for  interior  ballistics.  Calcula¬ 
tions  of  this  sort  are  capable  of  modeling  gun  performance  well  enough 
to  be  an  important  supplement  to  empirical  gun  tests  and  charge  firing 
trials.  In  addition,  some  understanding  of  gun  flash,  blast,  and  smoke 
have  been  gained  by  extensions  of  these  methods. 

We  are  greatly  indebted  to  the  personnel  of  the  U.S.  Naval  Weapons 
Laboratory,  Dahlgren,  Virginia  (principally  Dave  Bowen,  with  the  coop¬ 
eration  of  Frank  Kasdorf)  for  the  collection  and  interpretation  of  gun 
characteristics  and  firing  sequence  data  used  for  inputs  to  and  compari¬ 
son  with  our  calculations. 

Helpful  discussions  with  Alan  Shapiro  here  at  Rand  are  also  ac¬ 
knowledged  as  is  the  prior  interest  in  this  subject  by  Jack  Craig.  We 
would  also  like  to  thank  Wakichi  Asano  for  his  valuable  programing 
assistance. 

This  work  represents  a  progress  report,  since  the  authors  can 
foresee  further  Investigations  of  charge  characteristics,  energy  losses 
in  the  barrel,  and  other  non-ideal  effects,  as  well  as  a  more  detailed 
hydrody namical  model  of  the  smoke  and  flash.  Expressions  of  interest 
in  such  extensions  and  criticisms  of  the  work  presented  here  are  wel¬ 
comed  . 
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SUMMARY 


This  paper  considers  the  gas  dynamics  of  various  large  naval 
guns,  including  the  interior  ballistics  (charge  burning,  propellant 
gas  expansion  and  shell  acceleration)  and  the  nature  o£  the  gas  expan¬ 
sion  and  air  shock  beyond  the  muzzle.  Solutions  are  accomplished  bv 
a  numerical  program  which  integrates  the  partial  differential  equations 
appropriate  to  the  gas  dynamics  in  one  space  dimension.  The  equations 
of  state  for  the  gas  and  the  air  are  treated  in  some  detail,  as  is  the 
nature  of  the  charge  burn,  the  barrel  friction,  and  the  shell  inertia. 
The  dynamics  of  the  partially  burned  charge  gases  which  sometimes 
flash  burn  after  ejection  from  the  muzzle  are  treated  within  t lie  limits 
of  expansion  into  various  conical  geometries.  Comparisons  are  made 
with  observations  of  both  the  interior  ballistics  and  of  the  smoke, 
flash  and  blast  from  the  exhaust. 


Utr^ri  v*.ir  -■ 
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I_.  INTRODUCTION 

Tills  paper  introduces  an  accurate  numerical  method  for  obtaining 
the  hydrodynamic  properties  associated  with  the  propellant  gas  as  it 
moves  down  a  gun  barrel  and  exhausts  to  the  air.  Previous  analytical 
solutions  to  this  problem  have  found  it  necersary  to  make  a  number  of 
simplifying  assumptions  and  approximations,  some  of  which  limi:  the 
accuracy  and  range  of  applicability  of  the  results.  Even  the  very 
recent  treatment  of  high  speed  guns  by  A.  E.  Seigal^  takes  no  account 
of  the  nature  of  the  propellant  and  simply  assumes  that  it  burns  in¬ 
stantly  into  an  ideal  gas  with  constant  specific  heat.  Further,  Seigal 
assumes  that  the  shell  moveB  without  any  frictional  or  heat  losses. 

By  far  the  most  extensive  treatment  of  interior  ballistics  is  by 
J.  Corner.  ‘  He  suamarizefi  all  important  work  ir.  the  field  before 
1950.  The  "isothermal"  solution  is  the  simplest  and  it  assumes  that 
the  propellant  burns  according  to  a  specified  law  into  a  gas  which  ex¬ 
pands  isothermaliy  during  burning  with  a  constant  specific  heat  ratio 
which  is  increased  from  that  of  an  ideal  gas  to  account  for  neat  losses. 
The  shell  mass  is  increased  to  account  for  frictional  res'stanco. 

Generally,  in  Corner's  theory,  in  order  to  arrive  at  the  correct 
peak  breech  pressure,  muzzle  velocity,  pressure-time  history,  and  shell 
ejection  time  for  a  given  gun,  the  solution  is  made  to  depend  on  three 
propel  lsnt  parameters,  the  burning  law,  the  total  energy  released,  ar.ii 
the  rate  of  burning.  Corner  describes  a  more  detailed  solution  which 
takes  account  of  the  temperature  dependence  of  the  gas  and  includes  a 
number  of  other  small  corrections.  It  leads  to  very  complicated  ana¬ 
lytical  expressions  which  depend  basically  on  a  propellant  burning  law 
parameter  and  a  burning  iat<  parameter.  Corner  shows  a  few  cases  where 
good  agreement  with  experimental  -»’ues  can  be  obtained  by  proper  ad¬ 
justment  of  the  parameters.  Each  gun  ;st  be  treated  separately  and 
the  solution  is  very  sensitive  tc  the  various  ui justed  parameters- 

The  treatme:  vt  present  here  takes  good  account  of  many  of  the 
factors  in  the  problem  which  previously  suffered  from  c'^lifyinjr 
assumptions.  Using  a  full  partial  differential  formulation  ot  ib. 
dynamics  simplifies  sore  'actors  which  have  only  a  miner  effect  ni  lb- 


rwm  mi  bun* 


results  and  reduces  interior  ballistic  problems  to  a  single  free  param¬ 
eter,  namely  the  propellant  burning  rate.  Tills  parameter  is  obtained 
for  one  gun  by  fitting  to  experimental  data  and  the  solution  for  peak 
chamber  pressure,  muzzle  velocity,  pressure-time  history,  and  shell 
ejection  time  is  then  shown  to  be  accurate  to  the  order  of  ten  percent 
for  several  different  guns. 

Previous  descriptions  of  gun  flash  and  smoke  have  almost  entirely 
been  qualitative  in  nature  and  no  rigorous  mathematical  treatment  is 
knovm  to  exist.  The  hydrodynamical  solution  which  we  present,  while 
limited  in  Its  thoroughness  by  a  number  of  simplifying  assumptions, 
di.es  reproduce  quantitatively  a  number  of  properties  observed  in  actual 
guns,  such  as,  the  pressure  as  a  function  of  time  and  distance  from 
the  muzzle  and  the  shock  reheating  of  the  gas  outside  of  the  gun. 
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II.  INTERIOR  BALLISTICS 


In  thn  result:,  presented  'ere  are  examples  for  six  diflerent 

naval  Runs  ranging  in  si2e  from  5"/ 3b  caliber  to  lfa"/3()  caliber.  The 

/  >\ 

physical  characteristics  of  these  guns  are  listed  in  Table  The 

table  includes  such  parameters  as  chamber  volume,  propellant  weight, 
shell  weight,  and  shell  travel  distance.  All  the  propellants  used  in 
tnese  guns  are  essentially  identical.  They  consist  of  roughly  90/, 
nitrocellulose,  2-3/.  butyl  stearate,  1/.  ethyl  centralite  ,  14  lead 
carbonate,  17.  potassium  sulphate,  34  volatiles,  and  1-27.  moisture. 

They  come  in  the  form  of  cylindrical  grain  roughly  of  diameter  d 

o 

and  length  *  =  2dQ,  roughly  six  wel»  perforations  of  diameter 

Values  for  these  quantities  are  also  given  ir.  Table  1. 

In  this  analysis  of  interior  ballistics,  we  have  made  the  follow¬ 
ing  simplifying  assumptions:  (1)  each  gun  is  unworn  so  that  no  gas 
leaks  between  the  shell  and  the  barrel;  (2)  the  only  losses  in  the 
system  arc  frictional  in  nature;  (3)  the  pressure  drop  needed  to  over¬ 
come  the  skin  friction"  of  the  barrel  in  moving  the  gas  down  the 
barrel  is  negligible;  ( ~»)  gun  recoil  and  the  energy  expended  in  shell 
rotation  are  negligible;  (3)  the  chamber  can  be  replaced  by  j  cylin¬ 
drical  continuation  of  die  barrel  with  the  same  tocal  volume;  and  (6) 
the  propellant  grain  distritmees  itself  uniformly  hvtwecn  toe  breech 
and  shell  during  its  burning.  Tliese  assumptions  can  oe  shown  to  bn 
consistent  with  the  general  accuracy  of  our  method. ^  Other  conditions 
necessary  in  formulating  the  solution  are  discu-sed  below. 

MATHEMAT I CAL  FORMULATION 

The  equations  which  describe  the  interior  ballistics  part  of  the 
calculation  are  on  'dimensional  since  the  cross  section  is  assumed  to 
be  uniform  along  the  length  of  the  gun.  The  equation  cf  motion  in 
Lagrangian  form  is 

*  &  ill 

-s  the  velocity  of  tl»c  gas  in  meters; mi lli- 


-  si_ 

^  a 


where  u  u(x.t)  *  PAfc**! 

ot 
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second;  X(x,t)  is  the  Eulerian  coordinate  which  gives  the  position, 
at  time  t,  of  a  gas  element  that  was  initially  at  position  x,  the 
Lagrangian  coordinate;  m  ■  m(x,t)  =  /  p(x,t)dx  is  the  mass  of  the 

2  0  3 

gas  in  megagrams/meter“;  p  =  p(x,t)  is  the  density  in  megagrams/meter  ; 

3  16 

P  =  P(x,t)  is  the  pressure  in  jerks/meter  (1  jerk  =  10  ergs  or 

3  4 

1  jerk/meter  =  10  bars)  ,  and 


Q  =  Q(x,t) 


C^Am)2 

V 


3V 

3t 


2 


CjAm 

~v 


av  ,  ,  av  n 
H  I  for  sT  ‘  0 


(2) 


ana 


•  0  £or  f 


>  0, 


is  an  artificial  viscosity  pressure  introduced  to  remove  mathe¬ 

matical  discontinuities  which  would  otherwise  occur  whenever  compres¬ 
sions  lead  to  shock  waves  in  the  gas.  Typical  values  appropriate  to 

this  problem  are  C-  =  6.0  and  C,  =  0.5.  4m  is  a  mass  element  and 
i  /  3 

V  =  V(x,t)  =  l/p(x,t)  is  the  specific  volume  of  the  gas  in  meters  / 

megagram. 

The  equation  of  energy  conservation  is 


§§  -  -  (P  +  0)  ||  +  B  (P,f)  (3) 

where  E  =  E(x,t)  is  the  internal  energy  density  of  the  gas  in  jerks/ 
roegagram  and  B(P,f)  is  a  source  term  which  represents  the  rate  at 
which  energy  is  being  added  to  the  gas  by  the  burning  of  the  propellant. 
B  can  be  written  as  a  function  of  pressure  (P)  and  of  the  fraction  of 
propellant  already  burned  (f).^ 

For  me:  '-  propellants  the  burn  rate  in  energy  per  unit  time  (Bt> 
is  known  to  be  directly  proportional  to  the  exposed  propellant  surface 
area  (S)  and  the  gas  pressure  on  the  propellant  (P)  or 

Bt  =  cSP  (4) 

where  c  is  a  constant.  This  relationship  has  been  verified  for  nitro- 
cellulose,  '  although  at  pressures  near  atmospheric  P  should  be 
replaced  by  a  +  bp,  where  a  and  b  are  constants. 


For  the  cylindrical  grain  described  earlier,  neglecting  the  burn¬ 
ing  of  the  ends  and  the  webbing  as  a  first  approximation,  we  can  write 
S  as  a  function  of  f  and  other  known  constants.  For  the  fraction  of 
total  propellant  burned  as  a  function  of  time  we  can  write 

f  =  f(t)  =  1  -  ™d  -j/4-  =  1  -  (5) 

TTd0  i/4  dD 

Now  we  assume  that  the  propellant  of  total  weight  W  consists  of  n 

prop 

grains  each  of  weight  W  and  density  p  ,  so 

n  n 


W  =  nW  *  np  *'  =  n  p  ndn  S/4 

prop  n  Hn  n  Kn  0 


4  W  d 

S  =  nS  =  nndi  =  - ^2£_ 


so,  using  Eq.  (5)  , 


4W  d  4W 

B  ,  c  _ES!E- P  .  c  _E£1!  (1.f)l/2  p 

C  n  PnaO 


Thus  we  can  write  B(P,f)  in  energy  per  unit  mass  per  unit  time  as 


B(P, f )  -  s  (l-f)1/2  P  (9) 

4W  4 

where  s  *  c  — P--?-?-  —  =  ,01  m  /Mg-ms/d0.  The  numerical  value  .01  has 
p  d« 

n  (8) 

been  estimated,  independent  of  W  and  p^,  from  burning  rate  datav 

and  then  modified  to  give  a  better  fit  to  the  experimental  data,  as  will 

be  discussed  later.  Also,  we  assume  our  nitrocellulose  propellant  has 

a  calorific  value,  E  s  1000  cal/gm  -  4.2  jerks/Mg. ^ 

prop 

The  equation  of  motion  of  the  shell  is  obtained  by  using  the  ve¬ 
locity  and  pressure  of  the  gas  at  the  shell  boundary  as  it  moves  down 


the  barrel: 


du(«>h,t) 

sh  dt 


P(x9h.t)  A  -  Ffr 


where  A  is  the  cross  sectional  are^  of  the  barrel  and  F^  assumes  that 
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all  the  energy  lost  In  the  motion  of  the  shell  Is  caused  by  a  friction¬ 
al  force  acting  between  the  shell  and  the  barrel  which  is  directly 
proportional  to  the  pressure  on  the  shell  (appropriate  for  a  fluid  on 
a  deforming  plastic  seal  or  metal  shell) : 


Ffr  =  fr  ‘  P(xsh*t)  A‘ 


Gq.  (10)  can  be  rewritten  as  follows: 


du  PA  . ,  ,  .  PA 


where  M  ,,  =  8/rv  =  M  +  M  ,  . 

eff  (1-fr)  sh  sh  1-fr 


The  fraction  of  kinetic  energy  of  the  shell  going  into  frictional 
energy  loss  is  thus  fr/(l-fr).  Comer  suggests  that  this  fraction  is 
roughly  twenty  percent. 

Other  models  for  the  dissipative  force  were  tried  but  seemed  to 
be  much  poorer  in  terms  of  physical  soundness  and  resultant  solutions. 

One  such  alternate  scheme  assumed  the  frictional  force,  Eq.  (11),  was 
a  constant  and  another  subtracted  energy  from  Eq.  (3)  as  a  function  of 
the  kinetic  energy  of  the  shell. 

In  an  effort  to  derive  an  appropriate  caloric  equation  of  state 
for  nitrocellulose  we  have  relied  on  the  similarities  obvious  in  the 
Chapman -Jouguet  adiabatics  for  other  explosives ,  as  we  were  unable 
to  find  any  experimental  detonation  data  for  nitrocellulose.  In  Fig¬ 
ure  1  we  have  plotted  T-  1  vs  p  for  TNT  with  loading  density  pQ  - 
1.30  gm/cm  and  for  Pentolite  with  pG  «  1.65  gm/cm3.  By  fitting 


a  smooth  curve  to  this  data  in  the  region  of  interest,  naaciy,  p  £ 
.5  gm/cm  ,  we  obtain 


7(P)-1 


.09  +  2r 


.3  +  p4 


Nitrocellulose  propellant  generally  has  a  loading  de.-.Ly  c.  .7 
gm/cm  .  This  lover  initial  density  prumpted  a  somewhat  arbitrary 


I 


Density,  S  I  Mg  m3> 


Fig.  1— Caloric  equations  of  state 


modification  to  the  fit: 


r<P>- 1  - 


.09  + 


.3  + 


using  the  fact  that  the  value  of  y(p)-l  in  the  high  density  limit 
3 

(p»l  gm/cm  )  has  been  shown  to  be  linearly  related  to  loading  density. 
The  formulation  of  this  result  by  Miller,^1  based  on  explosives,  shows 
that  we  can  write 


r(»)- 1  -  3  +  g  Po/pc 


where  pQ  is  the  loading  density  of  the  explosive,  pc  is  its  crystalline 
density,  and  G  is  a  constant  determined  by  the  Morse  atomic  parameters 
associated  with  each  of  the  elements  composing  the  explosive.  An 
independent  determination  by  Fickett  and  Wood using  various  thermo¬ 
dynamic  relationships,  along  with  experimental  data  on  the  thermal 
equation  of  state  and  detonation  velocity  for  various  explosives,  leads 
to  empirical  results  for  y(°°)-l  similar  to  Eq.  (15).  The  transition 
from  the  high  density  region  to  the  low  density  region  of  the  gas  is 
obtained  by  using  Eq.  (14). 

Thus  our  caloric  equation  of  state  is 

2 

P  *  Ep[y(p)-1)  -  Ep  (-— -  -  A)  (16) 

.3  +  p 

For  low  densities,  p«.l  pn/cm^ ,  Eq.  (16)  describes  an  ideal  gas  of 
y*  1.30.  The  effect  on  the  results  of  varying  the  density  dependence 
of  y(P)*1  h**  investigated  and  is  reviewed  in  the  Results  and 

Discussion  Section. 

The  thermodynamic  description  of  the  gas  is  completed  by  using  a 
thermal  equation  of  state  of  the  form 


T  =  (H/R)  (P/p)  /  (1  +  /3(p) ) 


(17) 


where  R  is  the  Rydberg  gas  constant,  M  is  the  average  molecular  weight 
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o£  the  propellant  gas  (about  30  gm/mole)  ,  and  /3(p)  Is  a  "covolume" 

correction  ^to  the  ideal  gas  law,  which  is  roughly  linear  in 

3 

density  and  about  1.-2.  at  p  =  .5  gm/cm  ,  depending  on  the  exact  form 
of  the  correction  used. 

NUMERICAL  SOLUTION 

A  solution  for  the  hydrodynamic  functions  in  this  problem  results 
from  numerical  Integration  of  the  two  equations  of  motion,  Eqs.  (1) 
and  (12),  the  energy  equation,  Eq.  (3),  and  the  caloric  equation  of 
state,  Eq.  (16).  These  can  be  written  in  difference  equation  form  as 
follows: 


,n-l/2 


nf  1/2  n-1/2  r  n  n  . 

J  Uj  4m  LFj+l/2  >1/2  +  Qj+l/2 


nf  1/2  n-1/2  JL  p"  ) 

J  max  UJ  au  Hsh  Pj  max  U  tr} 


n-l/21  (18) 

Vl/2J 


(19) 


E)-l/2  EJ-l/2  l2(Pj-l, 


*  [2(Pj-l/2  +  Pj-l/2)  +  Qj-l/2]  (Vjtl/2  -  ^.l/2) 

(20> 


with 


nfl  _  „nfl  nfl  r  ,  nfl  .2-,/r  nfl  .2, 

P J - 1/ 2  E j -1/2  p J - 1/ 2  ^*09  +  j-1/2^  J/^*3  +  (pj-l/2) 

nfl  n  nfl/2  _  nfl/2 

J  j  °j  At 

nf  1/2  ^  2Cl(^l-l/2)  (V^l/2  ‘  V^-l/2) 

J'1/2  (V^l/2  *  ^-l/2><Atnfl/2)2 


(21) 


2C2  *1-1/2  I  <l/2'V?-l/2 

^1  +vn  77^172 
j-1/2  +  ”j-l/2 


(v1 


(23) 
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for  V 


n+1 


<  V. 


J— 1/2  J— 1/2 ' 


with  V 


n+1 

J-l/2 


n+1  Yn+1 

A  ~  J-i 

"“j-l/2 


(24) 


where  n  _  1  and  1  1  j  t  j  mas.  The  mass  elements  Am^  ^  are  fixed 
initially  and  remain  unchanged  throughout  the  solution.  Notice  that 
in  writing  the  source  term  of  Eq.  (20)  we  have  used  the  pressure  and 
energy  density  of  the  breech.  The  pressure  and  energy  density  are  not 
actually  uniform  throughout  the  gas,  but  as  discussed  later,  these 
approximations  have  a  negligible  effect  on  the  results. 

Unit  increments  in  the  superscript  n  represent  step-vise  incre¬ 
ments  in  time  and  unit  increments  in  the  subscript  j  represent  finite 
increments  in  the  mass,  corresponding  to  a  spatial  zone,  where  J  ■  1 

is  the  zone  at  the  beginning  of  the  chamber  and  j  ■  j  is  the  zone 

max 

adjacent  to  the  shell.  The  integration  is  performed  by  advancing  in 

time  (n)  with  the  step  size  determined  by  stability  conditions  which 

depend  on  the  local  sound  speed  in  the  gas  and  the  artificial  vlscos- 
(18) 

ity .  Then  for  each  step  in  time,  we  obtain  a  complete  hydrodyna- 

mical  description  of  the  propellant  gas  by  advancing  in  mass  zone  from 

J  «  1  to  J  ■  J  .  These  equations  yield  a  self-consistent  solution 
BAX 

for  u,  E,  p,  P  as  functions  of  time  and  position,  with  T  determined 
independently  from  p  and  P. 

By  using  the  loading  density  and  an  initial  energy  density  of 
the  propellant  assuming  atmospheric  temperature,  we  specify  gas  condi¬ 
tions  in  the  chamber.  The  constants  fr,  which  determines  the  effect¬ 
ive  mass  of  the  shell,  E  ,  the  energy /mass  released  from  the  pro¬ 
prop 

pellant,  and  s,  the  rate  of  propellant  burning  must  also  be  specified. 

From  previous  discussion,  we  know  fr  *  .2,  E  «  4.2  Jerks/Mg,  and 

prop 

s  *  .01/d..  These  values  are  used  initially  and  then  modified  to  ob¬ 
tain  the  best  possible  results.  The  effect  of  varying  these  parameters 
la  discussed  in  the  results. 


RESULTS  AND  DISCUSSION 

In  Table  2  we  present  for  each  of  the  six  guns  which  we  have 

studied  the  best  fit  values  for  the  three  quantities,  u  ,  the  shell 

muz 
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muzzle  velocity,  P^,  the  peak  breech  pressure  and  t^,  the  shell 

ejection  time,  along  with  the  corresponding  experimental  values.  In 

Figure  2,  we  present  calculated  and  experimental  pressure-time  curves 

for  the  5"/54  gun.  All  solutions  were  obtained  using  fr  -  0,  a.-*d 

E  “4.2.  We  used  s  *  1.48  for  the  5"/54  gun  and  then  computed 
prop 

s1  3  .01/d  as  an  estimate  for  the  other  guns.  The  actual  value  of 

°  i 

s  used  along  with  s  are  shown  in  Table  2.  The  results  show  agree¬ 
ment  between  the  best  fit  values  and  experimental  values  to  be  berter 
than  5*  in  most  cases  and  better  than  10T  in  all  cases. 

In  Table  3  we  present  a  summary  of  the  efficiency  of  guns.  Ti  rre 
are  two  standard  measures  of  efficiency.  The  first  is  called  "piezo¬ 
metric  efficiency"  and  is  defined  as  the  average  pressure  necessary  to 
accelerate  the  shell  to  the  observed  muzzle  velocity  divided  by  the 
peak  breech  pressure.  We  compute  the  average  pressure  P  as  follows 
from  Eq.  (13),  with  fr  -  0, 


P  A 


M  ,  a 
sh 


M 


,u  \ 

(  muz  ) 


sh  2  R 


tot 


or 


1/2 

sh  muz 


A  R 


(25) 


(26) 


tot 


where  R  is  the  combined  chamber  length,  L  ,  ,  and  shell  travel  length, 
tot  ch 

L  .  Piezometric  efficiencies  are  usually  about  .4  -  .5.  We  see 

trav 

from  Table  3  that  the  efficiencies  for  these  guns  are  somewhat  lower 
than  this. 

The  second  measure  is  "ballistic  efficiency,"  and  is  defined  as 

the  ratio  of  the  kinetic  energy  of  the  shell  at  muzzle  velocity  to 

the  energy  of  the  propellant.  We  shall  calculate  this  efficiency  in 

two  ways:  first  we  shall  use  the  full  energy/mass  of  the  propellant, 

E  *  4.2  Jerks/Mg,  and  second  we  shall  use  the  energv/mass  actually 
prop 

burned  at  the  time  of  si.wt  ejection,  E  (T  ,).  Both  these  efficiency 

prop  e) 

values  along  with  E  ^  (t  ,),  are  shown  in  Table  3.  The  tvnical  hal- 
prop  ei 

listic  efficiency  of  1/3^  '  lies  roughly  between  our  two  set*  o!  cal¬ 
culated  values.  The  maximum  possible  ballistic  efficiency  for  these 

guns,  obtained  if  the  charge  is  completely  burned  before  the  shot  moves 

{ Is*) 

appreciably,  can  be  shown  to  be  roughly  1/2. 


Breech  pressure,  P  lUilobors) 


I 


Fig.  2— Pressure-time  history  for  5"/ 54  gun 
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In  order  to  interpret  the  significance  of  these  results,  we  must 

discuss  the  changes  caused  by  varying  the  parameters  in  the  problem. 

These  changes  were  obtained  by  keeping  fixed  all  parameters  but  the  one 

whose  affect  we  are  studying,  and  then  obtaining  solutions  for  several 

values  of  this  parameter.  Changing  fr  from  .0  to  .25  decreases  u  ? 

by  less  than  4%  and  increases  P  by  about  25%,  with  the  changes  in 
1  max 

u  and  P  roughly  linear  with  the  changes  of  fr.  Changing  E 
muz  max  prop 

from  4.2  to  3.0  jerks/Mg  without  changing  s,  decreases  u  by  about 

3%  and  increases  P  about  1%.  One  reason  that  this  change  has  so 
max 

little  erfect  is  that  for  E  *  4.2  jerks/Mg  the  propellant  is  only 

about  60-75%  burned  at  i  . ,  and  lowering  the  value  of  E  to  3.0 

ej  prop 

Jerks/Mg  simply  increases  the  fraction  of  propellant  burned  at  t  .  . 

Lowering  E  still  further  would  then  begin  lowering  u  and  P 

prop  muz  max 

quite  strongly.  The  fact  that  the  nitrocellulose  propellant  we  have 

U3ed  is  only  partially  burned  at  indicates  that  it  most  likely  has 

a  calorific  value  less  than  1000  cal/gm.  Modifying  the  form  of  yG  )-l 

from  Eq.  (14)  so  that  it  goes  from  the  minimum  of  .3  at  low  density  to 

a  maximum  which  is  1.2,  increases  u  by  about  10%  and  P  bv  about 

muz  max  • 

20%.  The  shape  of  pressure-time  history  is  somewhat  less  sensitive  to 

changes  in  y(p)~1  but  if  the  maximum  value  is  raised  significantly 

above  1.0,  the  peak  in  the  curve  occurs  sooner  than  1/3  and  if  the 

maximum  value  is  much  lower  than  1.0,  then  ’•he  peak  occurs  at  roughly 

1/2  t^.  The  results  are  fairly  dependent  on  y(p)-l  but  we  have  chosen 

an  expression  which  we  believe  to  be  most  accurate. 

After  we  have  established  values  for  fr.  E  ,  and  v(p)-l,  as  dis- 

P 

cussed  earlier,  the  quantity  s  is  our  only  variable  parameter.  We  study 

the  effect  its  variation  has  on  u  and  P  for  each  of  the  guns. 

muz  max 

Using  the  fixed  conditions  fr  =  9.  and  E  =  4.2,  with  y(p)  given  by 

prop 

Eq.  (14)  and  with  B(P,f)  given  by  Eq.  (9),  a  curve  of  versus  s  is 

shown  in  Figure  3  and  a  cur^e  of  P  versus  s  is  shown  in  Figure  4. 

Combining  Figure  3  and  Figure  4  and  plotting  peak  breech  pressure 
versus  muzzle  velocity  in  Figure  5  it  is  seen  that  as  the  burn  rate 


increases  the  pressure  rises  much  more  rapidly  than  the  velocity.  This 
reflects  the  fact  that  as  the  burn  rate  increases  the  pressure  builds 
up  faster  to  a  larger  peak  value,  while  the  amount  of  energy  put  into 
the  gas  approaches  the  total  energy  in  ihe  charge. 
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Fig. 3— Muzzle  velocity  versus  burning  rate 
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Fig. 4— Peak  pressure  versus  burning  rate 
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Fig.  5— Peak  pressure  versus  muzzle  velocity 
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We  have  also  Investigated  using  burning  laws  cf  the  form  B(P,f)  * 
s  Pa(l-fY)  where  .8  <  a  <  .1  and  2  -  y  1  10  and  have  found  that  u 

muz 

and  P  are  not  particularly  sensitive  to  a  and  y  as  long  as  s  is  ad¬ 
justed  to  keep  the  burning  rate  roughly  the  same.  As  mentioned  earlier, 
we  have  also  made  calculations  using  individual  zone  pressure  and  energy 
instead  of  breech  pressure  and  energy,  and  these  yield  the  same  results 
as  long  as  s  is  adjusted  slightly  to  keep  the  overall  burning  rate  the 
same.  In  addition  we  have  investigated  burning  the  propellant  accord¬ 
ing  to  Eq.  (10)  only  in  the  chamber  part  of  the  gun.  Ve  have  found  that 
for  a  fixed  value  of  s, burning  only  in  the  chamber  vs.  burning  all  the 
way  between  the  breech  and  the  shell  leaves  Pmay  approximately  the  same 

and  decreases  u  by  about  102. 
muz 

As  a  further  study  on  the  nature  of  the  burning,  we  have  initially 
ignited  only  the  first  20*  of  the  propellant.  This  is  done  by  setting 
the  initial  energy  density  of  all  but  the  first  202  equal  to  zero  and 
then  letting  the  hydrodynamical  calculations  propagate  the  ignited  gas 
down  through  the  remaining  propellant.  The  result  of  this  is  to  create 
shock  waves  in  the  gas  which  propagate  between  the  shell  and  the  breech. 
The  magnitude  and  time  dependence  of  the  resultant  oscillations  in  breech 
pressure  depends  critically  on  the  fraction  of  propellant  initially  ig¬ 
nited  and  also  on  the  burning  law  for  the  remaining  propellant.  We  be¬ 
lieve  that  this  initial  burning  process  thus  accounts  for  the  oscilla¬ 
tions  which  are  observed  in  the  actual  pressure-time  histories.  However, 
since  we  do  not  know  the  precise  experimental  details  of  the  ignition 
process,  we  have  retained  the  original  assumption  of  uniform  burning  as 
a  good  first-order  approximation,  consistent  with  accuracy  of  our  other 
assumptions. 

As  a  final  study  we  have  investigated  the  effects  of  mass  ablation 
and  turbulant  heat  transfer  on  our  results.  These  energy  loss  effects 
are  treated  independently  of  the  frictional  losses  and  require  modifi¬ 
cation  of  the  equations  of  motion  and  energy.  In  the  Appendix  we  give 
a  detailed  description  of  the  new  equations  and  show  that  the  results 
are  only  slightly  affected  by  the  ch.-nges. 
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III.  FLASH  AND  SMOKE 

We  now  treat  the  propellant  gas  as  it  exits  from  the  muzzle.  The 

gas  composition  is  such  that  it  usually  forms  a  combustible  mixture 

with  the  surrounding  air  and  a  momentary  flash  of  bright  light  results 

from  a  shock- ini dated  burn.  Based  on  our  survey  of  available  experi- 

(3)  (20) 

mental  data  and  on  previous  studies,  the  flash  phenomenon  can 
be  divided  into  basically  three  region*.  The  first  two  regions  are 
muzzle  glow,  a  small  region  of  low  luminosity  at  the  muzzle,  and  the 
primary  flash,  a  high- luminosity  region  which  is  separated  from  the 
muzzle  glow  by  a  dark  region.  These  regions  are  usually  discernible 
only  for  small  guns.  For  the  large  guns  which  we  are  considering,  the 
third  region,  or  secondary  flash,  is  a  highly  luminous  continuation  of 
the  primary  flash  and  occupies  most  of  the  gas  volume. 

The  cause  of  the  flash  phenomenon  has  been  explained  qualitatively. 
Basically,  the  gas  expands  and  cools  very  rapidly  upon  leaving  the 
muzzle,  and  consequently  the  luminosity  diminishes  and  the  dark  region 
results.  This  gas  is  then  recompressed  through  a  decelerating  shock, 
which  raises  the  temperature  tr  ok  up  to  roughly  the  muzzle  temperature 
and  forms  the  luminous  primary  flash  which  is  relatively  small  and 
usually  lasts  for  a  few  milliseconds.  In  the  meantime  air  has  mixed 
with  the  unburned  hydrogen  and  carbon  monoxide  in  the  muzzle  gas,  form¬ 
ing  a  combustible  mixture.  The  exact  composition  of  the  muzzle  gas  is 

highly  dependent  on  such  parameters  as  the  initial  loading  density  of 

( 121 

the  propellant  and  the  final  density  of  the  gas,  besides  being 
dependent  cn  the  original  propellant  composition.  If  the  temperature 
of  the  recompressed  gas  is  above  the  ignition  temperature  of  the  com¬ 
bustible  mixture  or  if  a  piece  of  burning  propellant  penetrates  the  gcs 
cloud,  then  it  will  bum  and  fora  the  large  secondary  flash.  The  occur¬ 
rence  of  flash  is  often  Inhibited  by  the  presence  of  certain  substances, 
such  as  potassium  sulphate,  which  act  as  suppressants.  If  the  unburned 
propellant  gas  is  only  partially  ignited  or  not  ignited  at  all,  then 
it  appears  as  smoke  outside  of  the  gun. 

In  order  to  form  a  hydrodynamical  model  of  the  flash  phenomenon 
we  will  make  a  number  of  assumptions  which  simplify  the  calculations: 
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(1)  Che  muzzle  gas  emerges  Into  a  radially  symmetric  conical  section, 
which  Is  Initially  a  full  hemisphere  and  then  steadily  decreases  to  a 
very  narrow  angle  cone;  (2)  the  emerging  muzzle  gas  displaces  the 
surrounding  air  and  no  explicit  account  is  taken  of  any  mixing  between 
the  two;  (3)  the  unburned  propellant  remains  uniformly  mixed  throughout 
the  gas,  but  the  chemical  composition  of  the  gas  itself  is  not  deter¬ 
mined  or  used  in  the  calculations;  (4)  the  emerging  propellant  gas 
obeys  Eq.  (14)  and  the  surrounding  air  is  treated  as  an  ideal  diatomic 
gas  of  specific  heat  ratio  y  *  1.4.  Although  these  assumptions  are 
quite  restrictive,  we  still  expect  to  gain  certain  quantitative  informa¬ 
tion  about  the  nature  of  the  flash  phenomenon. 
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MATHEMATICAL  FORMULATION 

The  numerical  solution  now  separates  into  two  parts.  The  interior 
ballistics  are  still  governed  by  the  previously  described  equations  of 
motion  and  caloric  equation  of  state.  For  the  expansion  into  air  we 
have  the  equation  of  motion 

£  *  - r2  %  <f  +V  <27> 


where  r,  the  radial  distance  from  the  muzzle,  now  replaces  X,  the 
linear  distance  used  in  Eq.  (1)  ,  Mp  =  47rF/rp(r ,t)r2dr  =  P  r^  is 
the  total  mass  in  an  Fth  of  a  sphere  where  F  <  1,  nip  =  Mp/47r  =  ^  p 
r^  is  corresponding  mass  per  steradian.  We  can  write  F  =  sin2  ^ 
where  0  is  the  angle  between  the  axis  and  the  boundary  defining  the 
cone  into  which  the  gas  is  expanding,  and  Qp  is  the  modified  artificial 
viscosity 


<)r  -  QfU.O  ■ 


c  3<S>' 


where  and  are  constants  and  Amp  is  a  mass  element. 

The  energy  equation  outside  the  gun  becomes 

§  ‘  •  (p  +  V  5T  +  SR  (29) 

where  SR  is  a  source  term  which  represents  the  burning  of  the  gas. 
We  shall  arbitrarily  assume  that  SR  can  be  written  as  a  function  of 
gas  temperature  T  in  a  form  similar  to  Eq.  (9)  , 

SR  -  c(l  -  f)  1/2r  (30) 


where  c  is  a  constant  which  will  be  adjusted  so  that  the  gas  will  be 
completely  burned  about  10  milliseconds  after  first  exiting  from  the 
muzzle.  T  is  obtained  by  using  the  thermal  equation  of  state,  E<|.  (17), 
which  reduces  to  the  ideal  gas  law  at  the  low  densities  outside  ot 
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the  gun.  The  air  surrounding  the  gas  docs  not  burn  and  so  we  set 
SR  «  0.  in  that  region. 

The  caloric  aquation  o£  state  for  the  gas  is  given  by  Eq.  (16) 
and  the  equation  for  air  is 

P  -  -4Ep  (31) 

The  approximation  of  air  as  an  ideal  gas  is  valid  for  the  relatively 
low  temperatures  and  densities  which  occur  in  the  flash  region. 


NUMERICAL  SOLUTION 


We  can  now  write  Eq.  (27)  to  (31)  in  difference  equation  fora 

uf1'2  ■  ■r*'2  •  Km  -  - 1/2  *  oP;;;22  •  3  <»> 
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where,  in  the  gas  region. 
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where  Dq  Is  the  gui.  barrel  dlaaeter. 

The  solution  is  obtained  by  defining  the  origin  as  the  center  cf 
the  nuzzle  and  supplying  as  Initial  data  the  velocity,  pressure,  den¬ 
sity,  and  temperature  of  the  propellant  gas  as  a  function  of  negative 

distance  from  the  muzzle,  r"  =  x"  -  R  _  for  0.  *  X?  1  R  and  1  1  1 

j  J  tot  *  j  tot  J 

1  j  at  the  time  when  the  shell  reaches  the  muzzle.  The  numerical 

max 

integration  uses  Eqs.  (18),  (20),  (21),  (22),  (23),  and  (24)  for  the 
plane  geometry  region  r^  >  D0/,,  where  Eqs.  (38)  and  (39)  give  the 
appropriate  geometrical  factors  for  the  transition  between  the  two 
regions. 

The  initial  conditions  for  the  air  outside  the  gun  are  obtained 

4 

by  assuming  standard  temperature  (T  •  .03  x  10  *K)  and  pressure 

<f>  *  9.5  x  1<T5  jerke/m3). 

These  calculations  were  carried  out  for  the  5”/54  gun  as  a  repre 
sentatlve  case  and  one  for  which  we  have  comparable  experimental  data 
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RESULTS  AND  DISCUSSION 

A  survey  of  several  dozen  firings  of  the  5"/ 54  gun,  as  recorded 

(3) 

on  loom  film  taken  at  96  frames/sec,  shows  flash  and  smoke  character¬ 
istics  which  are  in  good  agreement  wit1-  the  general  description  given 
earlier.  The  gas  initially  emerges  into  a  full  hemisphere  (9  =  rr/2,  F 
3  1/2),  after  about  10  milliseconds  it  emerges  into  a  cone  of  angle 
0stt/3  (F *  1/4)  ,  and  afcer  about  20  milliseconds  it  emerges  into  a 
cone  of  angle  (F  ■  1/8)  ,  After  about  5°  milliseconds  the  gas 

is  emerging  straight  forward,  which  it  continues  to  do  until  it  is 
completely  exhausted  about  100  milliseconds  later. 

In  the  firings  where  there  is  no  flash,  the  initially  emitted 
gas  forms  a  large  cloud  of  smoke  roughly  two  meters  long  and  two  meters 
wide  which  drifts  away  from  the  muzzle  slowly,  and  through  this  is 
shot  the  gas  which  emerges  at  a  smaller  angle  and  forms  a  second  cloud 
of  smoke  further  dewnscream.  In  the  firings  where  there  is  a  flash 
it  is  possible  to  see  the  dark  region  of  expanding  gas  just  outside 
of  the  muzzle,  and  then  about  1. -2.  meters  downstream  the  shock  foms 
which  initiates  tise  flash  burn  within  the  first  1C  milliseconds.  The 
extent  to  which  the  flash  spreads  throughout  the  rest  of  the  gas  de¬ 
pends  quite  a  bit  on  the  amount  of  flash  suppressant  ir.  the  gas.  Some- 
c^*es  flash  is  only  a  small  region  of  luminosity  which  extinguishes 
Itself  very  soon,  and  sometimes  the  Hash  consumes  all  of  the  gas  and 
is  very  bright,  extending  over  several  meters. 

We  have  obtained  by  hydrodynamical  solutions  for  the  gas  and  air 
outside  of  the  gun  for  the  cates  F-  1/2  (U  -  -/ 2),  F  *  1/4  (0  =  -/3)  , 

*  »d  F  *  1/8  (0--/ 4).  To  show  t!*  format ian  of  the  shock,  we  have 
plotted  gas  temperature  versus  radius  at  a  fixed  time  for  each  of  the 
three  values  of  F  (see  Figure  6).  Note  !ov  the  gas  cools  very  rapidly 
as  it  comes  cut  of  the  muzzle  and  then  at  a  distance  of  1.5  -  2.  meters 
the  temperature  raises  back  up  to  its  initial  value  as  a  compression 
shock  is  formed.  The  unbumed  gas  is  then  ignited  and  flash  results. 
Note  that  the  location  of  the  shock  region  docs  not  change  great iy 
with  a  change  in  F.  The  fact  that  the  shock  M stance  at  a  given  time 
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is  greater  for  smaller  F  shows  that  the  gas  moves  faster  when  confined 
to  a  narrower  cone.  In  addition  to  temperature,  our  calculations  give 
the  gas  pressure,  density,  velocity  and  internal  energy  as  a  function 
of  position  and  time. 

In  Figure  7  we  have  plotted  the  gas  pressure  (assuming  F  -  1/2) 

at  a  fixed  distance  from  the  muzzle  (r  =  2.m)  versus  time.  Note  the 

sharp  rise  in  pressure  as  the  gas  first  arrives  at  2.m  and  the  gradual 

decrease  back  to  atmospheric  pressure.  This  pressure-time  history  is 

(3> 

in  good  agreement  with  corresponding  experimental  data,  'except  the 
experimental  pressure  curve  usually  drops  somewhat  below  atmospheric 
and  then  forms  a  second  peak  of  roughly  the  same  magnitude  as  the 
first  about  10  to  15  milliseconds  later.  Based  on  our  experimental 
data,  this  second  peak  is  due  to  the  flash  burning  of  the  gas  and  it 
does  not  occur  if  there  is  no  flash.  The  fact  that  we  do  not  reproduce 
a  second  peak  is  probably  due  to  the  fact  we  do  not  burn  the  gas  in 
the  air  properly,  which  would  require  knowledge  of  the  exact  chemical 
composition  the  gas  and  its  mixing  process  with  air.  Also  there  are 
small  instabilities  in  our  numerical  integration  procedure,  caused  by 
the  transition  from  plane  geometry  to  conical  geometry,  which  became 
larger  with  increasing  time.  These  cause  spurious  oscillations  in  the 
pressure  profile  which  would  most  likely  obscure  a  second  pressure  peak 
even  if  it  was  caused  by  our  burning  procedure. 

In  Figure  S  we  have  plotted  peak  pressure  versus  radius  using  our 

calculations  for  the  three  angles  9  =  tt/2  (F  =  1/2)  ,  0  =  tt/3  (F  =  1/4)  , 

and  9  =  tt/4  (F  =  1/8),  along  with  the  actual  experimental  data  measured 
(3) 

at  these  angles.  Note  that  we  have  made  the  assumption  that  the 
pressure  at  a  given  angle  is  obtained  by  using  the  value  of  F  corres¬ 
ponding  to  that  angle  in  the  numerical  integration.  We  feel  that  this 
is  a  rough  approximation  to  the  actual  pressure  distribution.  Although 
the  calculated  pressures  are  greater  than  the  corresponding  experimental 
values  by  about  50%,  they  all  lie  on  straight  lines  with  almost  the 
same  slope,  thus  correctly  reproducing  the  adiabatic  expansion  of  the 
gas  as  it  moves  away  from  the  muzzle.  The  fact  that  we  have  had  to 
rely  on  separate  calculations  using  different  conical  angles  in  order 


ig.  7— Pressure-time  profile  for  5"/34  gun  at  r  =  2m 
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to  simulate  the  angular  dependence  of  the  pressure  points  up  a  major 
limitation  in  our  one -dimensional  conical  calculations. 

A  final  fact  which  comes  out  of  our  calculations  is  >:hat  the  time 
for  the  pressure  peak  to  travel  from  a  radius  r  =  2m  to  r  *  5m  is 
roughly  6ms,  which  is  in  good  agreement  with  the  experimentally  meas¬ 
ured  time. 
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IV.  CONCLUSIONS 

We  have  obtained  a  good  quantitative  model  for  the  interior  bal¬ 
listics  of  large  naval  guns.  Our  results  are  limited  by  a  number  of 
simplifying  assumptions  but  they  are  valid  within  these  limitations. 

Our  major  Improvement  over  previous  treatments  is  the  use  of  a  more 
realistic  equation  of  state  for  the  propellant  gas.  We  have  shown  how 
sensitive  the  interior  ballistics  are  to  a  number  of  parameters.  With 
more  detailed  experimental  data  we  improve  these  parameters  somewhat 
and  also  investigate  the  replacement  of  the  various  ideal  assumptions 
by  the  actual  physical  process. 

The  flash  and  smoke  are  much  more  difficult  to  model  in  detail. 
Although  our  treatment  accounts  fairly  well  for  the  pressure  distribu¬ 
tion  and  shock  reignition  of  the  gas,  we  need  a  much  better  knowledge 
of  the  mixing  of  the  gas  and  air  and  the  flash  burning  mechanism.  A 
number  or  the  mathematical  features  of  the  hydrodynamics  need  to  be 
Improved.  The  use  of  two-dimensional  calculations  is  necessary  to  pro¬ 
vide  for  the  angular  dependence  of  the  pressure.  Alternate  methods  of 
treating  the  transition  region  between  plane  and  conical  geometry  should 
be  studied  in  order  to  improve  the  stability  and  accuracy  of  the  calcu¬ 
lations  in  that  region. 

We  have  established  the  validity  of  our  calculations  for  large 
naval  guns,  but  in  order  to  show  their  general  validity  they  should 
be  applied  to  a  variety  of  other  guns. 
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APPENDIX 


A  study  of  some  gun  barrel  effects  which  have  a  small  but  finite 
effect  on  the  gas  dynamics  can  be  made  by  generalizing  the  Lagrangian 
equations  of  motion  to  include  the  effects  of  heat  and  mass  transfer 
due  to  ablation  and  gas  turbulence  processes  during  the  interior  bal¬ 
listics. 

The  equations  of  motion  are  recast  in  an  "almost-Lagrangian"  form 

which  allows  each  of  the  zones  to  increase  in  mass  as  wall  ablation 

(21) 

occurs.  Following  the  derivations  of  Crowley,  we  obtain  an  equation 
of  continuity 
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where  m  is  the  mass  flux  (mass/area/time)  and  S  is  the  surface  area 
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per  unit  mass.  The  momentum  equation  becomes 
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where  uy  is  the  velocity  of  the  mass  entering  a  zone  and  x^  is  the 
shearing  stress  at  the  wall.  The  energy  equation  becomes 

~  -  -  (P  +  Q)  ~  +  a  Sm  [<u-u  )2/2  +  E  -  E]  (A-3) 
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where  the  second  term  on  the  right  hand  side  gives  a  rate  of  energy 
addition  associated  with  the  mass  addition  and  E^  is  the  specific  in¬ 
ternal  energy  of  the  mass  entering  a  zone.  The  third  tern  comes  from 
work  done  on  a  mass  element  by  wall  surface  forces,  and  the  fourth 
term,  H,  is  the  rate  at  which  heat  is  supplied  to  an  element  by  turbu¬ 
lent  heat  transfer. 

The  following  expression  for  the  mas^  flux  is  used: 


„  C  (yPu)/(y-l) 
q  h _ 

*  "  E  +  i (PV+E) 
q  v 


Sru 


n  +  E 


v/ (PV+E) 


(A-4) 


where  is  the  dimensionless  coefficient  of  heat  transfer  (the  Stanton 
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number),  .Is  Che  sum  of  the  enthalpy  per  unit  mass  required  to  heat 
the  wall  material  to  vaporization  temperature,  the  specific  heat  of 
fusion  and  the  specific  heat  of  vaporization,  and  n  is  the  turbulent 
transpiration  coefficient.  The  source  term,  representing  the  heat  lost 
from  a  zone  due  to  the  turbulent  ablation  process,  is  vritten  as 

H  -  -  q  S  -  -  Cuo(PV+E)uS  (A-5) 

On  IB 

The  shearing  stress  can  be  written  as 

Tw  ‘  1  Cf  2  u2  *  CH  “  “2 

where  C,  is  the  dimensionless  coefficient  of  skin  friction  and  C  ■  ^  C, 
f  (2n  H  2  f 

by  the  Reynolds  analogy. 

To  obtain  a  numerical  solution  we  now  put  Eqs.  (A-l)  to  (A-6)  in 
difference  equation  form.  Using 
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and  substituting  into  Eqs.  (A-2)  ar.d  (A-3)  we  obtain 
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where  we  have  set  u  ,  *  0. 

v  ( 22 3 

Ue  have  followed  Crowley  and  Glenn  '  in  evaluating  the  constants 
for  a  steel  shell  moving  down  a  steel  barrel: 

CH  -  .002,  o  -  .2,  Ey  -  10  jerks/Mg,  Ey  -  12.5  jerks/Hg. 

To  test  the  sensitivity  of  our  results  to  these  values  we  have  varied 

ami  Ey,  which  have  the  most  effect  on  the  results. 

Results  show  the  maximum  breech  pressure,  P  ,  lowered  less  than 

n&x 

21  in  all  cases  and  the  muzzle  velocity,  u  ,  lowered  by  IX  for 

muz 

-  .002,  *  10  jerks/Mg;  21  for  -  .004,  E^  »  10  jerks/Mg; 

10%  for  »  .02,  E^  *  10  jerks/Mg  and  4%  for  *  .002,  ■  20  Jerks/Mg 

and  raised  bv  2%  for  C  *  .002  and  E  ■  5  jerks/Mg.  We  thus  hav-  -hown 
that  the  size  of  the  ablation  and  turbulent  heat  transfer  corrections 
is  small. 
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